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Abstract 

The ground state energy and pairing gap of the interacting Fermi gases calculated by the ah initio 
stochastic method are compared with those estimated from the Bardeen-Cooper-Schrieffer pairing 
Hamiltonian. We discuss the ingredients of this Hamiltonian in various regimes of interaction 
strength. In the weakly interacting {l/akp << 0) regime the BCS Hamiltonian should describe 
Landau quasi-particle energies and interactions, on the other hand in the strongly pairing regime, 
that is l/akp ^ 0, it becomes part of the bare Hamiltonian. However, the bare BCS Hamiltonian 
is not adequate for describing atomic gases in the regime of weak to moderate interaction strength 
— oo < l/akp < such as akp ~ —1. 
PACS: 05.30.Fk, 03.75.Ss, 21.65.+f 
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The superfluid state of alkali fermion gases such as ^Li, is analogous to the super- 
conducting state found in the electronic systems 111 and it has been studied theoretically 

Q n 

[2| and experimentally p||. Recent experimental progress includes the detection of the Bose 
Einstein Condensate (BEC) state of the bound Fermi pairs J, 0,0, El as an evidence of the 
predicted BCS-BEC crossover. 

The interaction among ^Li atoms can be attractive and requires preparation of the atoms 
in different internal quantum states |1) and |2). The Feshbach resonance 0,0,^3 between 
the atoms in these states can modify the characteristics of the collision. We assume that the 
partial densities of the atomic species are the same and the temperature is low (T ^ 0). The 
interactions can be characterized by the s-wave scattering length a which can be tuned by 
the externally applied magnetic field. The Hamiltonian of the interacting atomic fermions 
can be written in the configuration space as 

^2 N 

^'"^'^ = -|rEv? + Enr.O (1) 

p=l 

where the index i is for species |1) particles while i' is for species |2) particles. In the 
momentum space, the Hamiltonian can be written as 

^ are ^ ^ a^,k'^«T,k + 9 E %,k',qCtl,k+qCt2,k'-q'^2,k'ai,k (2) 

a,k ^ k,k',q 

with a\ 1^ and Oo-.k being particle operators. This is the so-called bare Hamiltonian. 

On the other hand, in terms of the quasi-particles, we have the so-called Landau Hamil- 
tonian of the form 



-^Landau ^ + ^ _ ^t^^^^^^ + ^ ^ <k',q4,k+q4,k'-qC2,k'Cl,k , (3) 

(T,k k,k',q 

where and Co-,k are quasi-particle operators. En = energy of the 'normal' ground state. 



= quasi-particle energy spectrum, and fic = Ep. In general, we have ^^k,k',q ^k,k',q- 

Further simplification yields BCS pairing Hamiltonian that can also be written in two 

ways 

'^BCS = E «f7,k«fT,k + E %,k''^l,k«2 -k«2 -k'fll.k' (4) 

o-,k k,k' 
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and 

n^^^^t"" = i^Tv + E [^k - /ic] <kc.,k + E <k'4,k4,-kC2,-k'Ci,k' (5) 

cr,k k,k' 

The BCS approach is to restrict the interaction to the time-reversed pair (k,— k) of different 
species. 

Ground states of Hamiltonian (Eq have been obtained using the stochastic method 



known as Green's Function Monte Carlo (GFMC) 



111 ll^ where variational degrees 



of freedom were introduced to deal with the fermion sign problem. Their energies are 
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121, a 



considered as close upper bounds to those of the exact ground state. In the Ref 
finite short range cosh-function potential rather than (5-function potential was used. The 
range of the potential is ~ |ro where tq is the unit radius (Ivrrgp = 1). From the results of 
the lowest order cluster calculations known as LOCV Q] (Lowest Order Constrained 
Variational) it appears that this finite range potential is a good approximation for the zero 
range potential in the l/akp < regime (see FigHJ. 

The original BCS technique was to use 7i^^^"" (EqEl) and solve it variationally using 
^Bcs- The solution will give the superfluid energy Es = + Ep, where the pairing energy 
Ep = —Nq^. In general, it is difficult to obtain analytically Ti^^g"'^ from 7i^^^. However, 
for weak potential strength, that is 1/akp < —1, we can map the 'H^p^fs by 
making the substitution 

(6) 

+ const. (7) 

(8) 

In this case we can get Es = E^ + Ep = Epenz + Ep (see Fig El). This projection of the 
bare Hamiltonian into the Landau Hamiltonian is useful only in the weak interaction limit 
in which Epenz — Epc has most of the interaction effect and \Ep\ « \Epenz — Epc] as seen 
in the FigH 
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On the other hand, Leggett 
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ig solved the H^'^s (Eq ^ with the condition that 



the density remains constant with the chemical potential /ic adjusted accordingly. This 
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TABLE I: Comparison of Leggett results with n'^^cs vs GFMC. The unit of energy is EpG 
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We notice Ep = 1.67EpG- We define AE = Escs-Leggett — Egfmc- We notice that while there 
is considerable discrepancy in the energies, the gaps are in reasonable match for l/akp < —1/3. 
Errors are in the last digit except for Agfmc where the relative error ~ 10%. 



method can be applied in all regimes of interaction. The interaction of the Hamiltonian 



assumes zero range v^y., 



g = = volume) adequate for the dilute regime where 

the potential range R and a are much less than tq ~ as well as the 'intermediate' regime 

we can draw the normalization 



where R « tq « a. From = |fk| = ■ 

condition. Going to the continuum limit and expressing in the units of Ep — , 
a set of two equations 



1 



dee'/' 



1 



+ A2 



A2 



4 
' 3 

71 

akf 



— , we have 



(9) 



(10) 



where Eq comes from subtracting the equation for the scattering length a (see Ref |l 
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FIG. 1: Comparison of E/N calculated using different methods. Both finite range 'cosh' potential 
and 5-function like potentials are considered for the LOCV calculations. 



from the gap equation 



ink dk 



1 



(12) 



^(e(fc)-/.,)2 + A2 

The fic and A are solved simultaneously. The solutions are given in the Table HI In this 
table, the energy per particle EBcs-ieggett = E/N was estimated using /Xc and A, and the 
expression 

(13) 



E/N = —J2'2eM^ - ^u^Vk 
k 

with the usual definitions of A = — (^I^Mk^^k and l^ikP = 1 — |^k| 

k 



In the Fig |21 we compare the normal phase low density Lenz expansion 
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(14) 
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FIG. 2: Comparison of E/N and A calculated using BCS-Leggett equations and the stochas- 
tic GFMC method. The EpG is subtracted from all energies. Second order El^uz and Ep = 
-Nq ^bcs -Landau ^^^^ ^g^^ ^^ot Eg = E Lenz + Ep. Es has good match with Egfmc up to 
l/akp ~ —1. 

of first and second order of akp with BCS-Leggett and GFMC results both from the 
bare potential. The Lenz expansion is considered exact in the low density regime for 
the interacting Fermi gas in normal phase. From Fig |2l it is obvious that the expansion 
diverges for 1/akp > —1 . Es = Ei^nz + Ep has good match with the GFMC results in the 
regime l/akp ^ —1. Here the effect of pairing is small in Egfmc (the difference between 
the Slater node and BCS node solutions lies within the statistical uncertainties) thus 
\Ep\ « lEienz — Epcl is a reasonable assumption although A is clearly non zero. In this 
regime we notice that Epcs-Leggett is distinguishably higher than Egfmc and energies of the 
normal phase low density expansion. At akp = —1, AE = Epcs-Leggett — Egfmc = 0.19EpG 
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(see Table H)). This is a consequence of having the anomalous density u^Vk small and 
fic ~ Ep. In fact, when u^Vk ~ the "^bcs ~^ ^fg- However, we can see that the 
usual Hartree-Fock term for the normal phase is missing (Eq IT^ in the energy expres- 
sion from the pairing Hamiltonian. Thus we have Escs-Leggett much higher even than Eienz- 

We interpret this as a consequence of using pairing Hamiltonian instead of the full bare 
interaction Hamiltonian. The pairing Hamiltonian becomes a poor model for the atomic 
gas in the interacting regime with — oo < 1/akp < 0, in particular around the moderate 
interaction strength akp ~ —1- This is in sharp contrast to the context in which the 
original Landau-BCS formalism was introduced that was the weak coupling approximation 
in a broad range -oo < 1/akp < -1. Egpmc, EBcs-Leggett, and E^enz + Ep converge in 
the afci? ~ regime, that is the trivial free Fermi gas limit, where the Hartree-Fock term 
becomes effectively zero. 

On the other hand, the time reversed pairing (k,-k) assumption becomes less relevant 
once the interaction is strong enough for the particles to form loosely bound pairs in the sea 
of many fermions. This can be seen as Ai? becomes smaller (= {).lbEpG) in the strongly 
interacting regime 1/akp = and 1/akp = +0.1. In the 1/akp > 0.3 region, EBcs-Leggett 
and Egpmc apparently reverse back to the diverging behavior (Fig^). But as shown in the 
comparison of the LOCV energies, the range of the model potential becomes inadequate 
to approximate the 5-potential as the size of the bound pairs become ^ ro and ~ R. We 
argue that the actual Egpmc with short range potential would lie closer to EBcs-Leggett 
than the current finite range calculation shows. The bound fermions condense in the k = 
state. Thus U^Zm (EqEJ ^ ^bcs (EqlH) and results of two models should match. 

As for the pairing gap, both BCS-Leggett and GFMC results seem to be in rea- 
sonable agreement in the whole 1/akp < region considering that statistical errors of 
^GFMC ~ 10%. The reasonable match of A while a poorer match for E /N is not surprising 
given the fact that the chemical potential /ic is greatly modified in this region, /i^ goes 
from ~ Ep = = 1.67EpG at 1/akp ^ —1 to the /ic < for 1/akp ^ 0.5 where zero 
momentum excitation is possible and BEG is achieved. 
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In conclusion, we have tested the regimes of vahdity of the BCS pairing Hamiltonian 
in the study of fermion particles interacting with bare short-range two-body potential. 
We notice that the pairing assumption is generally not valid when bare potential is used 
in a broad range of the weakly interacting regime — oo < l/akp < 0, while the original 
quasi-particle BCS formalism was introduced to describe the superfiuid precisely in this 
region. We notice considerable discrepancy in the energy, however the gap is predicted 
with reasonable accuracy at akp ~ —1. Pairing correlation is less relevant in the trivial 
(free Fermi gas) and the tightly bound pair (l/akp > 0) limits. In fact, it can be shown 
that GFMC calculations with both Slater and pairing nodes converge to the same value 
(molecular energy per particle Emoi/'^) in the extreme of this limit. This work has been 
supported in part by the US National Science Foundation via grant PHY 00-98353 and 
PHY 03-55014. The authors thank useful comments from Prof. G. Baym of UIUC and J. 
Carlson of LANL. 
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